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ABSTRACT 


ant F-type star from NASA's Transiting Exo- 
Rj, and an orbital period of 7.245991 0.00092 
90.58 

-0.57 


1 
0.088 
—0.035 


My, resulting in a bulk density of 14.1*17 g cm ?. This makes it one of the few massive giant planets with extreme density and lies in 
the transition mass region of massive giant planets and low-mass brown dwarfs, an important addition to the population of less than 
five objects in this mass range. The eccentricity of 0.325 + 0.020 and an orbital separation of 0.0888 + 0.0010 AU from its host star 


and metal enrichment of the planet (Zp /Zstar) of aot. 


v] [astro-ph.EP] 21 Mar 2023 


—4 1. Introduction 


ft Massive giant planets (4 to 13M;) have always been debated 
whether these objects should be classified as planets or Brown 


— dwarfs (BDs)(Chabrier et al.|2014j {Spiegel et al.|2011j Schlauf- 


: [man|2018). There are a few indirect ways to discriminate be- 
e tween the massive giant planets and the low-mass BDs. One of 


them is based on the deuterium burning mass-limit, which states 
that an object should not be massive enough to sustain deuterium 
e » fusion at any point in its life to be classified as a planet. The up- 
per mass limit for this deuterium fusion was calculated to be 
«13M, for objects of solar metallicity (Boss et al.|2005), re- 
gardless of their formation channel. However, the objects with 
less than 13M; share quite common "nature" with 13M; objects, 
irrespective of what they have been called. That's why this def- 
inition based on a clear-cut mass limit between BDs and plan- 
ets has caused disagreements (Chabrier et al.[2014), and various 
suggestions have been made to reshape it. The study by [Spiegel] 


(2011) proposed that deuterium burning may vary from 11 
to 16M ;, depending on the object's helium and other metal con- 


tent. Some other studies have recommended increasing the upper 
mass limit to ~ 25M; based on the "driest" region of the brown 
dwarf desert (Pont et al. 005 Anderson et al. 011). 
In another prospect, (2015) provided a new det- 


inition and preferred the objects in 0.3-60M ; to be called Gi- 


. 


rXiv:23 


e 


suggest that the planet is likely undergoing high eccentricity tidal (HET) migration. We find a fraction of heavy elements of 0.13 


+0.05 
—0.06 


Detection of such systems will offer us to gain valuable insights into the 
governing mechanisms of massive planets and improve our understanding of their dominant formation and migration mechanisms. 


Key words. stars: individual: TOI-4603 — planetary systems — techniques: photometric — techniques: radial velocities 


ant Gaseous planets as they follow a particular sequence in the 
mass-density diagram of all the known planets, sub-stellar ob- 
jects, and stars (see figure1, Hatzes & Rauer](2015)). They do 
not see any abrupt changes in the mass-density diagram for the 
objects in 0.3-60 M; and suggested irrespective of any forma- 
tion scenario, these objects should fall under the same general 
class of objects, i.e., planets. However, recently IAU proposed a 
working definition of exoplanets(Lecavelier des Etangs & Lis-| 
which states that in addition to the 13M, mass limit, 
the system should have a mass ratio with the central object below 
the L4/Ls instability (M/M central < 2/05 + V621) ~ 1/25). 
Some researchers favor the basis of formation mechanisms 
to distinguish massive giant planets from BDs. Theoretically, 
two formation mechanisms dominate the literature: core accre- 


tion (Pollack et al.|1996| generally followed by low-mass gi- 
ant planets (Mp < 4M), |Schlaufman| (2018)) and Disk insta- 


bility 1997| generally favored by massive giant planets 
as well as low-mass BDs). However, the dominating mecha- 


nism for planet formation depends on the disk mass and host 
star metallicity conditions, i.e., their initial environmental con- 
ditions (Adibekyan|2019). Therefore, it is not clear how to trace 
the formation history of a planet from the current understanding, 
and this definition is also inadequate and problematic. Hence, the 
detailed characterization of more massive giant planets and low- 
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mass BDs will enhance our knowledge of the processes involved 
in planet formation and offer more insight into the transition re- 
gions of these objects. 

One frequently debated aspect about the close-in massive gi- 
ants or giant planets is whether they are formed at their present- 


day short orbits or migrated from further out orbits (Batygin & 
[Stevenson|2010] [Baruteau et al.|2014). The common belief is 


that these planets form beyond the ice line and then migrate in- 
ward through various mechanisms to their present location. Mi- 
gration of a planet to a close-in orbit occurs via torques from the 
proto-planetary disk (gas disk migration) or gravitational scat- 
tering due to another planet or star (HET migration). Eventually, 
due to tidal forces, its orbit is circularized and shrunk (see[Daw-] 
and references therein for a detailed re- 


view). Nevertheless, recent models suggest the in-situ formation 
of close-in giant planets is also feasible and 
show that the inner boundary of short-period giant planets and 
their period-mass distribution could be consistent with predic- 
tions for in situ formation (Bailey & Batygin{2018). That's why 
which of these three scenarios predominated is still up for dis- 
cussion, but a combination of these processes likely contributed 
to the current close-in giant planet population. 

In this letter, we report the discovery of TOI-4603 b, a close- 
in massive giant planet in the overlapping mass-region of BDs 
and planets. The subsequent sections discuss all the observa- 
tions, analysis, and results. 


2. Observations 
2.1. TESS Observations 


TESS observed the star TOI-4603 (HD 245134) in three sec- 
tors 43, 44, and 45. AII the observations were made with the 
two-minute cadence mode nearly continuously between Septem- 
ber 16, 2021, and December 02, 2021 (~ 74 days time-span) 
with a gap of ~ 5.5 days due to the data transferring from the 
spacecraft. Light curves were produced and analyzed for tran- 
sit signals by the Science Processing Operations Center (SPOC: 


2016), consisting of Simple Aperture Photome- 
try (SAP) and Pre-search Data Conditioning Simple Aperture 


Photomeiry (PDCSAP: 
fluxes. These light curves are publicly available at the Mikulski 
Archive for Space Telescopes (MAST[] The SPOC pipeline de- 
tected ten transits with a depth of ~ 1020 ppm, an orbital period 
of ~ 7.24 days, and a duration of ~ 2.04 hours. We adopt the 
median-normalized PDCSAP fluxes for further analysis that are 
additionally detrended by fitting a high-order polynomial over 
out-of-transit data using the lightkurve package 
[Collaboration et al.[2018). The normalized TESS light curve for 


TOI-4603 is shown in Figure[1] The Target Pixel Files (TPFs) of 
TOI-4603 generated with tpfplotter (Aller et al.|2020) for all 
B.6 


the observing sectors can be found in Figure 


2.2. High resolution imaging 


To assess the possible contamination of bound or unbound close 


companions on the derived planetary radii (Ciardi et al.|2015), 


we observed the TOI-4603 with near-infrared adaptive optics 
(AO) imaging at Palomar Observatories. The observations of 
TOI-4603 were made with the PHARO instrument (Hayward 


2001) behind the natural guide star AO system P3K 


"https://mast.stsci.edu/portal/Mashup/Clients/Mast/ 
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Fig. 1: The normalized PDCSAP light curve for TOI-4603 is plotted 
with respect to time in the upper panel, and time in days from mid- 
transit is plotted in the lower panel. The 2-min and 10-min binned data 
points are represented by blue and orange dots, respectively. The black 
line represents the best-fitted transit model using the EXOFASTv2 (see 


SectionB.3). 


(Dekany et al.|2013) on November 21, 2021, in a standard 5- 


point quincunx dither pattern with steps of 5" in the narrow-band 
Bry filter (A, = 2.1686; AA = 0.0326 um). Each dither position 
was observed three times, offset in position from each other by 
0.5" for a total of 15 frames, with an integration time of 5.665 
seconds per frame, respectively, for total on-source times of 85 
seconds. PHARO has a pixel scale of 0.025" per pixel for a total 
field of view of ~ 25". The AO data were processed and analyzed 
with a custom set of IDL tools. The science frames were flat- 
fielded and sky-subtracted and then combined into a single im- 
age using an intra-pixel interpolation that conserves flux, shifts 
the individual dithered frames by the appropriate fractional pix- 
els, and median-coadds the frames. The final resolutions of the 
combined dither were determined from the FWHM of the point 
spread functions: 0.117". The sensitivities of the final combined 
AO image were determined by injecting simulated sources az- 
imuthally around the primary target every 20? at separations of 
integer multiples of the central source's FVHM 
[2017). The brightness of each injected source was scaled until 
standard aperture photometry detected it with 5c significance. 
The resulting brightness of the injected sources relative to TOI- 
4603 set the contrast limits at that injection location. The final 5c 
limit at each separation was determined from the average of all 
of the determined limits at that separation, and the uncertainty 
on the limit was set by the RMS dispersion of the azimuthal 
slices at a given radial distance. The final sensitivity curve for 
the Palomar data is shown in (Figure B.1); no additional stellar 
companions were detected. 


Gaia Assessment 


The Gaia Renormalized Unit Weight Error (RUWE) is a metric 
similar to a reduced chi-square, where values that are x 1.4 in- 
dicate that the Gaia astrometric solution is consistent with the 
star being single, whereas RUWE values 2 1.4 may indicate an 
astrometric excess noise, possibly caused by the presence of an 


unseen companion (e.g., Ziegler et al.|2020). TOI-4603 has a 
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Gaia EDR3 RUWE value of 0.998, indicating that the astromet- 
ric fits are consistent with the single-star model. 


2.3. Spectroscopy 
2.3.1. Radial Velocities with PARAS 


RV observations were obtained using the PARAS spectrograph 
coupled with the 1.2m telescope at PRL Gurushikhar Observa- 
tory, Mount Abu, India. PARAS is a fiber-fed echelle spectro- 
graph with a resolving power of R=67000, a wavelength cover- 
age of 380—690 nm. A total of 27 spectra were acquired between 
January 11, 2022, and November 02, 2022, using the simulta- 
neous wavelength calibration mode with Uranium-Argon (UAr) 
hollow cathode lamp (HCL) as described in [Chakraborty et al.] 
and [Sharma & Chakraborty] (2021). The exposure time 
for all the spectra was 1800 s leading to SNR per pixel of ~ 9-18 
at the blaze wavelength of 550 nm. More details on observations 


and data analysis can be found in|Chakraborty et al.|(2014). The 


reported uncertainties are measured in the same way as described 


in |Chaturvedi et al.|(2016| 2018); [Khandelwal et al.|(2022). All 


the RVs and their respective errors are listed in Table 


2.3.2. Radial Velocities with TRES 


We obtained 13 observations between November 3, 2021, and 
September 16, 2022, using the Tillinghast Reflector Echelle 
Spectrograph (TRES; |Fürész|2008) on the 1.5m Tillinghast Re- 
flector telescope on Mount Hopkins, AZ, USA. TRES is a fiber- 
fed echelle spectrograph with a resolving power of R=44,000 
and operating in the wavelength range 390-910 nm. The spec- 
tra were obtained in a sequence of 3 observations surrounded by 
ThAr calibration spectra, and then the median combined to re- 
move cosmic rays. The average exposure time was 290 s result- 
ing in an average SNR per resolution element of 54.2. The spec- 
tra were extracted using procedures outlined in [Buchhave et al.| 
(2010), and multi-order relative velocities were derived by cross- 
correlating the strongest SNR observed spectrum order-by-order 
against all of the remaining spectra. RVs acquired with TRES 
spectra with their respective errors are listed in Table[A-2] 


3. Analysis 
3.1. Spectroscopic Parameters of TOI-4603 


We used the Stellar Parameter Classification (SPC; 
Buchhave ct al 2014) to de- 
rive stellar parameters from TRES spectra. SPC cross-correlates 
an observed spectrum against a grid of synthetic spectra based 
on Kurucz atmospheric models (Kurucz]1992). Using 12 of the 
13 spectra that passed the quality flag based on SNR, we derive 
Te =6243 + 50 K, log g, of 3.94 x 0.10 cgs, [m/H] of 0.22 x 0.08 
dex, and vsini of 25.70 +0.50 km s^. 

We also obtained high SNR spectra (70 per resolution ele- 
ment at 550 nm) of the 1200 s with Tautenburg coude echelle 
spectrograph (TCES) installed at the 2m Alfred Jensch tele- 
scope, Thüringer Landessternwarte Tautenburg, Germany. TCES 
is a slit spectrograph with a resolving power of R=67000 and a 
wavelength coverage of 470—740 nm. For details of the observa- 


tions, one can refer to |Guenther et al.|(2009). The spectra were 


extracted using the IRAF and used to compute stellar param- 


eters with zaspe package (Brahm et al.|2017). It yields Tog of 
6273: 101 K, log g, of 3.73+0.26 cgs, [Fe/H] of 0.34 £0.04 dex, 


and vsini of 23.18 € 0.37 km s^! through comparison against a 


grid of synthetic spectra generated from the ATLAS9 model at- 
mospheres (Castelli & Kurucz[2003). The stellar parameters ac- 
quired from TRES and TCES spectra are within the error bars. 

Our analysis shows that the TOI-4603 is a metal-rich, F- 
type sub-giant star. We also calculate the star's rotation period 
by computing Generalized Lomb-Scargle periodogram (GLS; 
on the out-of-transit TESS PD- 
CSAP light curves and find it to be 5.62 x 0.02 days which is 
comparable to the rotation period (assuming i=90) derived using 
vsini (section B.1) and stellar radii (section p.3). A less signifi- 
cant peak at ~2.28 days was also observed in the periodogram, 
which may be quasi-periodic and related or unrelated to half of 
the rotational period signals. Prewhitening the 5.62 days signal 
did not eliminate the 2.28 days signal, possibly suggesting it 
originated from another active region on the staller disc. How- 
ever, further analysis of the 2.28 days signal is beyond the scope 
of the current work. 

We also inspected the star for solar-like oscillations in the 
star. We first calculate the expected frequency of the maxi- 
mum oscillation amplitude (vmax) using the above calculated Teff 
and logg, with the seismic scaling relation (Lund et al.[2016), 
which yields Vg, «700 wHz. Since this value is smaller than 
the Nyquist frequency for the 2-min (~4166 uHz) cadence data, 
TESS photometric data is well-suited for identifying the oscilla- 
tions. We analyzed the oscillation signals using the 1lightkurve 
package and manually studied the power density spectra of the 
same TESS light curves but could not detect any significant 
solar-like oscillations. 


3.2. Periodogram Analysis 


Independent of photometry, we search for periodic signals in RV 
data from both spectrographs, PARAS and TRES, using the GLS 
periodogram. These RVs are corrected for the instrumental offset 
prior to analysis. The periodogram is shown in the panel 1 of Fig- 
ure[B-2] Here, we calculate the false alarm probability (FAP) of 
signals using equations given in [Zechmeister & Kürster| (2009) 
and find the most significant signal at 7.24 days (marked with 
vertical red line in Figure [B2]. This period is the same as es- 
timated from transit data (see Section ET}. The signal gives a 
FAP of 0.007% at 7.24 days using a bootstrap method over a 
narrow range centering this period, robustly confirming the pe- 
riodic signal in our RV data set. The other significant signals in 
the RV periodogram vanish after removing the 7.24-days peri- 
odic signal using a best-fit sinusoidal curve into the datasets (see 
in panel 2, residual periodogram). The spectral window function 
is shown in panel 3. As a diagnostics of stellar activity indica- 
tor and stellar contamination from nearby stars, we compute the 
periodogram of bisectors in panel 4 and find no statistically sig- 
nificant signal of stellar activity in the data sets. 


3.3. Global modeling 


We constrained the system parameters with simultaneous model- 
ing and fitting of the RVs from PARAS and TRES and the TESS 
light curves using the publicly available EXOFASTv2 
package. The software incorporates the differential 
Evolution Markov Chain Monte Carlo (MCMC) technique with 
the Bayesian approach to explore all the given parameter space. 

EXOFASTV2 uses a combination of spectral-energy distribution 


(SED; 2016) modeling, the stellar evolutionary 


models; generally MESA Isochrones and Stellar Tracks (MIST) 
isochrones (Choi et al.|/2016} [Dotter|2016); and the prior pa- 
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Fig. 2: The obtained RVs from PARAS and TRES are plotted with re- 
spect to ~7.24 days orbital phase. The best-fitted RV model with EXO- 
FASTv2 (see Section 3.3) is represented by the red line, and residuals 
between the best-fit model and the data are shown in the bottom panel. 


rameters to constrain the host star parameters. We performed 
the SED fitting for TOI-4603 using the broadband photome- 
try from Tycho BV ANTR, SDSSgri, APASS DR9 
BV(Henden et al.|2016), 2MASS JHK (Cutri et al.|2003), ALL- 
WISE W1, W2, W3, and W4 E ke in 
Table [A.T] We imposed Gaussian priors on the Teg and [Fe/H], 


determined from spectral analysis of the TCES spectra. Along 
with that, we also provided a Gaussian prior on parallax from 


GaiaDR3 (Gaia Collaboration et al.|2022) and enforced an up- 
per limit on the V-band extinction of 1.59 from |Schlafly & 


(2011) dust maps at the location of TOI-4603. The 
SED fitting uses Kurucz's stellar atmospheric models 


[1979), and the resulting best-fitted SED model with broad- 
band photometry fluxes is shown in Figure [B.3] Within the EX- 
OFASTv2, the MIST evolutionary tracks are used to provide 
better estimates for the host star parameters. The most likely 
MIST evolutionary track from EXOFASTVv2 provided the age of 
1.64*00 Gyr (see Figure B.4). The adopted stellar parameters 


are Tep =6264*33 K, logg,—3.810*0 07^ dex, [Fe/H]-0.342*0 07 
dex, M,=1.765 +0.061 Me, and R,—2.738*0055 Ro. All the pa- 
rameters are summarized in Table[A.3]|along with their 1-c- un- 
certainty. 

The simultaneous fitting of the RV and transit data is done by 
keeping all the parameters (like b, i, Rp, a K c, e) free and only 
providing starting values of P and T, given by the TESS QLP 
pipeline. The transit model of was used 
for light curve fitting while the RV data is modeled with a stan- 
dard non-circular Keplerian orbit. We used the default quadratic 
limb-darkening law for the TESS passband and the limb dark- 
ening coefficients (u; and u2) were calculated based on tables re- 
ported in |Claret & Bloemen|(201 1I) and [Claret| (2017). We used 
42 chains and 50000 steps for each MCMC fit that are further 
diagnosed for convergence using built-in Gelman-Rubin statis- 
tics [Ford|2006). The transit and RV 
data with their best-fitted models using EXOFASTV2 are plotted 
in Figure [1] Figure D] and Figure [B.5] In RV data, we also fit- 
ted a long-term RV trend (y) and found it to be —0.14 + 0.18 
ms day | (Table [A.3]. which may not significant due to its 
higher uncertainty. All the planetary parameters obtained by EX- 
OFASTV2 are reported in Table[A.3] 
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4. Results and Discussion 
4.1. TOI-4603 b in context 


We find the mass and radius of TOI-4603 b as 12.89*058 M; and 


-0.57 
1.0427) 078 Rj, respectively, transiting an F-type sub-giant star 


in a 7.24599*000077 days orbit. The discovery of TOI-4603 b is 
a substantial contribution as it is in the overlapping mass region 
(11 to 16M;; [Spiegel et al.[2011) of massive giant planets and 
low-mass brown dwarfs (BDs) based on the deuterium burning 
mass limit. As per the IAU definition, for solar metallicity, the 
deuterium burning mass limit is 13M, 
& Lissauer|2022). However, this limit depends on other factors, 
such as the abundance of helium and initial deuterium, and on the 
metallicity of the invoked model. For example, for three times 
the solar metallicity, 10% of initial deuterium can start burning 
at 11M; (Spiegel et al.|201 1). Assuming the metallicity of TOI- 
4603 b to be the same as that of its parent star, i.e., 0.342700 
dex, the companion here would have initiated deuterium fusion, 
loosing on its first criterion to be called a planet. However, ac- 
cording to the second criterion, TOI-4603 b has a mass ratio of 
0.007, with the host below the L4/L5 instability ( «1/25), which 
is in favor of it being called as an exoplanet. Finding the explicit 
nature of the astrophysical body in this mass range, whether it is 
a planet or a BD, can be an ambiguous task (see|Schneider et al.| 
(2011) for a detailed overview). For many in the field, includ- 
ing See aN and references therein), do not consider 
the deuterium burning mass limit as a strict boundary to distin- 
guish planets and BDs. There are other studies that suggest the 
upper mass limit for a planet, in particular, a gas-giant planet 
should be 25M, 
and, in some cases 60M ;(Hatzes & Rauer|2015). Since 
TOI-4603 b, according to most of these definitions, qualifies as 
a gas-giant, we would rather call it a planet in our current work. 

We, hereby, present the mass vs. density plot of transiting 
gas-giant planets and BDs that have mass and radius with a 
precision better than 2596 with reported mass ranges between 
0.25M; (lower mass limit for the gas giants from [Dawson &| & 
to 85M, («0.08Mo) in Figure [3] To date, there 
are a total of 5310 confirmed exoplanets, out of which 1569 
exoplanets' masses have been determined'| Here we focus on 
the transiting giant planets (0.25-13M j), which leave us with 
477 transiting giant planets, and of these, 35 are massive giant 
planets (Mp > 4M j)*| We plot the Mp=13M,, the deuterium fu- 
sion mass limit for solar metallicity, as a vertical dotted line and 
shaded area for planet-BD overlapping mass-region. As seen in 
the Figure, there have only been three such close-in (a«0.1 AU) 
transiting objects (HATS-70 b: and XO-3 b: 
Johns-Krull et al. (2008) discovered in this mass range, includ- 
ing our work. This makes TOI-4603 b an important addition in 
the context of all the known giant planets. 


4.2. Internal structure 


We estimate the heavy element content of TOI-4603 b using 


the method described in (2021). Given the prop- 


erties of TOI-4603 b, we estimate the planetary radius obtained 
with the evolution model completo21 
and compare it with the observed radius. We assume that all the 
heavy elements are homogeneously mixed in the envelope and 


'"http://exoplanet.eu/ 
*https://www.astro. keele.ac.uk/jkt/tepcat/ (2011) as of 


November 16, 2022 
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Fig. 3: Planetary density as a function of planetary mass for transiting 
giant planets and brown dwarfs (0.25-85 M). The shaded area repre- 
sents the overlapping mass region of massive giant planets and brown 
dwarfs based on the deuterium burning limit, and the dotted lines are at 
Mp=13M, and Mp=85M,, respectively. The position of TOI-4603 b is 
denoted by the magenta dot. 


are modeled as water with the Equation Of State (EOS) of water 


ANEOS ('Thompson|1990;|Mordasini|2020). Similarly to[Thorn-] 
gren & Fortney| (2018); |Komacek & Youdin|(2017), we do not 


include a central core. The envelope is coupled with a semi-gray 
atmospheric model, and hydrogen and helium (He) are modeled 
with SCvH EOS with a He mass fraction 
Y 20.27. We use a Bayesian framework to infer the internal lu- 
minosity of the planet, which matches the planet's radius given 
its mass and equilibrium temperature. The internal luminosity is 
governed by a linear uniform prior, and the content of heavy el- 


ements is informed by the|Thorngren et al.|(2016) relation. We 


find that the planetary radius is well reproduced with a fraction 
of heavy elements of 0.137005, As noted in (2021), 
the prior on the internal luminosity has an effect on the final in- 
ternal luminosity; however, the two values of heavy elements are 
compatible within 1o. From this fraction of heavy elements in 
the envelope, we can derive the metal_enrichment of the planet 


Zp|Zstar=4.2*5° (as done in $4.3 from|Ulmer-Moll et al.|(2022)) 


and the total mass of heavy elements M,-532*75 Mo. We in- 
clude in Appendix D the posterior distribution of the fitted pa- 
rameters. 

TOI-4603b is a scientifically interesting object for studying 
the processes of planet formation at the transition between mas- 
sive giant planets and BDs. proposed two 
populations of giant planets with masses above and below ~ 4M; 
in their study. Specifically, their finding suggests that the forma- 
tion of lower-mass giant planets may be related to core accre- 
tion and have metal-rich hosts. In contrast, higher-mass planets 
may form through disk instability mechanisms and orbit stars 
with lower average metallicity values. Moreover, 
established this theory by finding that planets with Mp 
< 4M ; preferentially orbit metal-rich hosts, unlike planets with 
Mp > 10M; do not have this trend. With the high metallicity 
([Fe/H]20.342*0-05? dex), TOI-4603 b does not follow this trend 


—0.040 
and does not support the existence of any breakout point at 4 Mj, 


as suggested by |Adibekyan| (2019). It demonstrates that regard- 


less of the metallicity of the host star, a massive giant planet can 


be formed via any process (Adibekyan|2019). 


4.3. Eccentricity of TOI-4603 b and tidal circularization 


The orbit of TOI-4603 b is found to be eccentric (e=0.325 + 
0.020). Different processes, such as secular interactions, planet- 
planet scattering, planet-disk interactions, and high-eccentricity 
tidal migration, explain the orbital evolution of giant planets 
(see Section 2 of for more details). 
We plot the observed population of the transiting giant planets 
(0.25M; < Mp < 13Mj) in eccentricity and semi-major axis pa- 
rameter space (similar to [Dong et al.] (2021)) in Figure [4] us- 
ing the TEPcat databasdil The shaded area indicates the region 
where planets could have undergone HET migration following 
the constant angular momentum tracks. The boundary of this re- 
gion is determined by the Roche limit and the tidal circulariza- 
tion timescale, respectively, and is defined as a=0.034—0.1 AU. 
The position of TOI-4603 b indicates that its orbit is undergoing 
HET migration. 

Based on the giant planets’ eccentricity distribution (Fig- 
ure (4). planets with orbital periods between 3 and 10 days have 
a wider range of eccentricities (0.2 < e < 0.6) than those with 
shorter periods (e < 0.2). HET migration is the most favor- 
able explanation for these moderate eccentricities, implying that 
these eccentric giant planets are in the process of tidal circu- 
larization. We also observe circular and eccentric giants at the 
same orbital periods, reason being, circular giant planets started 
their migration earlier than eccentric giant planets or have more 
efficient tidal dissipation effects. Some low eccentricities may 
be due to other formation channels like in situ formation or 
disk migration. Furthermore, mostly all eccentric giant plan- 
ets orbit metal-rich stars, whereas circular giant planets orbit 
both metal-poor and metal-rich stars (Figure ff}. Given the well- 
known correlation between the occurrence of giant planets and 
stellar metallicity, established 
that eccentric giant planets primarily orbit metal-rich stars. Their 
findings support HET migration via planet-planet gravitational 
interaction. Being a metallic host, and eccentric orbit of TOI- 
4603 b is consistent with this trend. Moreover, [Kervella et al.] 
found that the TOI-4603 has a widely separated (~1.8 
AU) BD companion (Mp « 20.52M j) in its orbit. This BD com- 
panion of the TOI-4603 system may provide an explanation for 
this eccentricity. We also calculated the shortest tidal circular- 
ization timescale (Teir) of 8.2 Gyr (for Q=10°;|Adams & Laugh-| 
[lin|2006), greater than the star's current age determined from 
this work. So, as per tidal evolutionary theory, the orbit of TOI- 
4603 b has not been circularized, which is consistent with our 
observations. 


5. Summary and Future prospects 


This work presents the discovery and characterization of a tran- 
siting giant planet around a subgiant star TOI-4603 at an or- 
bital period of 7.24599*009522 days and was initially identi- 
fied as an exoplanet candidate using transit observations by 
NASA’s TESS mission. Further, we complemented the TESS 
data with the ground-based observations with PARAS/PRL, 
TCES/TLS, TRES, and PHARO/Palomar instruments. Based 
on the global modeling of the TOI-4603 system, the host star 


is found to be metal-rich ({Fe/H]=0.342*)3) dex), subgiant 


(log g.=3.8107}035 g cm”), F-type (Te#=6264793 K) star that 
has a mass, radius, and age of 1.765+0.061 Mo, 27381002. 


Ro, and 1,547039 Gyr, respectively. The planet TOI-4603 b has 
a mass of 12.89*058 Mj, a radius of 1.042*0028 Rj, and ec- 


centricity of 0.325 + 0.020 with an equilibrium temperature of 
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Fig. 4: The orbital eccentricities of all the transiting giant planets (0.25M; < M < 13M) are plotted with respect to their semi-major axis (a) in 
AU. The datasets are taken from TEPcat database} and planets with a precision better than 25% in eccentricities are considered. The gray region 
represents the path of high-eccentricity migration with a range of 0.034-0.1 AU in the final semi-major axis and is set by the Roche limit and the 
circularization timescale, respectively. The giant planets are color-coded according to their host's metallicity. Triangles represent planets with P<3 
days, and diamonds represent planets with 3«P«10 days, and circles represent P> 10 days. The position of TOI-4603 b in the figure is marked 


with an arrow. 


1677 + 24 K. It is one of the most massive and densest transit- 
ing giant planets known to date and a valuable addition to the 
population of less than five massive close-in giant planets in the 
high-mass planet and low-mass brown dwarf overlapping region 
(11M) «Mp < 16M) that is further required for understanding 
the processes responsible for their formation. 

TOI-4603 is a rapid rotator (vsini=23.18 + 0.37 km s7!) 
and relatively bright star (V=9.2), well suited for Rossiter- 


McLaughlin effect (R-M; 1924} [McLaughlin] [1924] 


study and helpful in measuring the projected stellar obliquity of 


planets. The calculated RM semi-amplitude (Ohta et al.|2005) 


for the projected spin-orbit angle (4) between 0° and 90° is 
6.4 m s^! and 31 m s^! , respectively. The detection of the RM 
effect for TOI-4603 is possible by observing precise RVs us- 
ing moderate-sized telescopes (2.5-4 m aperture); for example, 
PARAS-2 (Chakraborty et al.[2018) at 2.5 m telescope, PRL is 
well suited for this work. 
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Table A.1: Basic Stellar Parameters for TOI-4603 


Parameter Description (unit) Value Source 
Œ J2000 Right Ascension 05:35:27.82 (1) 
Ô 2000 Declination +21:17:39.62 (1) 

Ho PM in R.A. (mas yr^!) 0.102 + 0.021 (1) 
Us PM in Dec (mas yr7!) -22.866 + 0.011 (1) 
m Parallax (mas) 4.4613 + 0.0195 (1) 
G Gaia G mag 9.0831 + 0.0027 (1) 
T TESS T mag 8.6554 + 0.0062 (2) 
Br Tycho B mag 9.964 + 0.026 (3) 
Vr Tycho V mag 9.273 x 0.019 (3) 
B APASS B-mag 9.915 + 0.03 (6) 
V APASS V-mag 9.42] +0.15 (6) 
g SDSSg mag 9.968 x 0.23 (6) 
r SDSSr mag 9.310 + 0.18 (6) 
i SDSSi mag 8.976 + 0.04 (6) 
J 2MASS J mag 8.089 + 0.020 (4) 
H 2MASS H mag 7.788 + 0.047 (4) 
Ks 2MASS Ks mag 7.786 + 0.017 (4) 
WI WISEI mag 7.718 + 0.028 (5) 
W2 WISE2 mag 7.744 + 0.02 (5) 
W3 WISE3 mag 7.761 + 0.02 (5) 
WA WISE4 mag 7.933 + 0.198 (5) 
Ly Luminosity (Lo) 9.74 [9.65, 9.80] (1) 
Tet Effective Temperature (K) 6189 [6185, 6193] (1) 
logg Surface gravity (cgs) 3.805 [3.801, 3.818] (1) 
[M/H] Metallicity (dex) -0.236 [-0.239, -0.232] (1) 
M, Mass (Mo) 1.752 + 0.088 (1) 
Ry Radius (Ro) 2.722 + 0.136 (1) 
Age Age (Gyr) 1.98 [1.73, 2.22] (1) 
Other Identifiers: 
HD 245134? 
TIC 437856897? 


TYC 1309-1102-1? 
2MASS J0535278242117396* 
GaiaEDR3 3402980516507429888! 


Note: The metallicity of TOI-4603 reported by Gaia is different 
from our spectroscopic analysis (see section 3.1). 
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Table A.2: Radial Velocity measurements of TOI-4603. 


BJDrpg Relative-RV o-RV BIS c-BIS  EXP-TIME Instrument 
Days ms! ms! ms! ms! S 
2459591.245018 1301.87 57.30 -2207.44 268.35 1800 PARAS 
2459592.214846 1155.71 62.52 1217.54 199.04 1800 PARAS 
2459619.190349 791.20 52.30 -92.82 249.29 1800 PARAS 
24596 19.224340 744.52 82.10 413.35 278.35 1800 PARAS 
2459619.269151 589.99 83.52 236.22 255.28 1800 PARAS 
2459647.154375 91.81 58.45 -22.06 128.41 1800 PARAS 
2459647.190552 -18.86 54.08 -331.87 91.57 1800 PARAS 
2459648.118842 584.50 84.86 -1495.12 285.70 1800 PARAS 
2459650.113942 1211.70 66.83 -1008.67 149.20 1800 PARAS 
2459650.191793 1082.77 86.56 -1048.42 206.49 1800 PARAS 
245965 1.119375 160.29 62.99 -451.66 147.03 1800 PARAS 
2459651.150425 251.80 58.78 -991.27 247.10 1800 PARAS 
2459673.160786 -35.14 82.03 259.99 224.52 1800 PARAS 
2459676.145032 -172.45 79.06 -2088.29 231.24 1800 PARAS 
2459678.118658 1374.26 99.56 -745.09 421.19 1800 PARAS 
245988 1.366954 1586.09 59.53 -1033.86 151.34 1800 PARAS 
2459881.390763 1655.28 70.41 -1160.88 206.22 1800 PARAS 
2459882.340429 831.30 5621 -414.69 100.21 1800 PARAS 
2459882.363486 838.82 52.98 248.93 439.60 1800 PARAS 
2459882.498600 390.78 50.84 -417.80 230.02 1800 PARAS 
2459883.322491 16.65 48.72 -1154.95 86.43 1800 PARAS 
2459883 .346278 -116.47 44.56 -393.03 170.19 1800 PARAS 
2459884.297863 -287.89 56.43 -2149.04 106.84 1800 PARAS 
2459884.321626 -328.58 48.01 -1984.25 99.19 1800 PARAS 
2459885 .323514 -226.52 96.60 -4272.63 284.06 1800 PARAS 
2459886.321280 25.85 55.50 -3578.19 330.80 1800 PARAS 
2450886.418752 1.46 68.72 2356.22 159.67 1800 PARAS 
2459521.904947 -509 130 - - 90 TRES 
2459525 .893170 1012 69 - - 180 TRES 
2459526.860655 1436 54 - - 450 TRES 
2459604.831871 234 99 - - 270 TRES 
24598 19.997097 -552 46 - - 360 TRES 
2459820.998886 -396 58 - - 180 TRES 
2459824.001461 1193 96 - - 195 TRES 
2459824.981289 -130 50 - - 400 TRES 
2459829 .013338 0.00 78 - - 720 TRES 
2459830.019859 1210 116 - - 300 TRES 
2459836.992174 470 59 - - 210 TRES 
2459837.976241 1474 95 - - 180 TRES 
2459839.014168 540 78 - - 240 TRES 
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Table A.3: Priors along with Median values and 68% confidence interval for TOI-4603 from EXOFASTv2. The N and U represent the Gaussian 


and the Uniform priors, respectively. 
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Parameter Units Adopted Priors Values 

Stellar Parameters: 
M..... Mass(Mo)........... cence eee ee - 1.765 x 0.061 
Rosas, CRAdiSRel: coii ue tret ors vts - 23380. 
Dolo Luminosity (o) etext lesen, = 10.40+0-65 
Ps... Density (CgS). sensere ap UE - 0.121 1706677 
logg.. Surface gravity (cgs)................ - 58100] 
Tef... Effective Temperature (K)........... N(6169, 128) 626475 y 
[Fe/H]  Metallicity (dex) ................00. - 0.34270 
Age... Age ( Gyi lsc Set Ries tetro elena - 6403 
EEP..  EqualEvolutionary Point ........... - 299.1555 
Ay....  V-band extinction (mag). ............ U(0, 1.5965) 0.225008 
CsEp. SED photometry error scaling ....... - 3,647085 
vsini.. Projected Rotational Velocity (km s^!) - 23.18 + 0.37 
(E Parallax (mas).............02e eee e ee N(4.4613, 0.01947) 4.462 + 0.020 
d... Distance (po). ........ssseeesesseess — 224.12 +0.99 
Vine ais RV slope (m/s/day) ................. — —0.14+ 0.18 

Planetary Parameters: b 
Pus. Period (days)... oec eth does - 7.24599 oop 
Rp... "Radius Op) sait couette irr eis - 1.042075 
Tc.... Time of conjunction (BJDrpp)....... - 2459549.1260 + 0.0014 
ü en Semi-major axis (AU)............... - 0.0888 + 0.0010 
ER Inclination (Degrees) ............... - 80.2123" 
Eis Eccentricity 22: ee Rs - 0.325 + 0.020 
w,.... Argument of Periastron (Degrees)... . - 20.4746 
Teq... Equilibrium temperature (K)......... — 1677 +24 
Mp... Mass (M)... aonana annaa - 12807275 
Kis RV semi-amplitude (m/s)............ E 96273 
logK.. Log of RV semi-amplitude .......... - 2.983 + 0.016 
Rp/R,. Radius of planet in stellar radii ...... - 0.03917 
a/R...  Semi-major axis in stellar radii ...... - 6.97 x 0.14 
eed Transit depth (fraction).............. - 0.001528+9:000001] 
Depth. Flux decrement at mid-transit ....... - 0:001598 ae 
Ti4... Total transit duration (days) ......... — 0.1189 + 0.0022 
b. eus Transit Impact parameter ........... - 0.952176 p049 
pP.... Density (cgs). .... 0... cece cece eee - 14.17 
loggp. Surface gravity ..............suuuu. - 4,4690036 
(F)... Incident Flux (10? erg s! cm^?)..... - 1.622803 
Tp.... Time of Periastron (BJDrpg)........ - 2459548.363*0075 
ECOS ^ nore S ON e er RU Ur SOR UE — 0.303 + 0.019 
QU E ——M m 0.113 x 0.027 
Mpsini Minimum mass (My) ............... - 12705927 
Mp/M,  Massratio ..........0. 000 c ee - 0.006988 38028 

Wavelength Parameters: TESS 
uj.... linear limb-darkening coeff ......... 0.237 + 0.050 
u2 .... quadratic limb-darkening coeff ...... 0.318 + 0.050 

Telescope Parameters: PARAS TRES 
Yret--. Relative RV Offset (m/s) ... ......... 2967 147*2 
Cogo IWJditer(m/s)i o seeccerir gs 95175 185*8* 
o3... RY Jitter Variance n.on 91007500 35000*10000 


Transit Parameters: TESS (TESS) 
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0?.... Added Variance ................sss. 0.0000000151 + 0.0000000018 
Fo.... Baseline flux ...................us. 1.0000094 + 0.0000024 
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Appendix B: Figures 
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Fig. B.1: Palomar NIR AO imaging and sensitivity curves for TOI-4603 
taken in the Bry filter. The images were taken in good seeing conditions, 
and we reach a contrast of ~ 7 magnitudes fainter than the host star 
within 0."5. Inset: Image of the central portion of the data, centered on 
the star. 


Period (days) 


* 9S S o o S S 
€ 9 x 9 D 


—— Residual RVs 


—— Window function 


0.8 
Frequency (cycles per day) 


0.0 02 04 06 1.0 


Fig. B.2: The GLS periodogram for the RVs, residual RVs, window 
function, and bisector slope of TOI-4603 is shown in panel 1, 2, 3, and 
4 (upper to lower) respectively. The primary peak is seen at a period 
~ 7.24 days (red vertical line), consistent with the orbital period of the 
planetary candidate obtained from photometry. The FAP levels (dashed 
lines) of 0.1% for all the periodograms are shown in the legends of the 
panel 1. 
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Fig. B.3: The spectral energy distribution of TOI-4603. 
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Fig. B.4: The most likely MIST evolutionary track from EXOFASTv2 
for TOI-4603 represnted by solid black line. The other two green dashed 
lines show the evolutionary track of 1.58 Mo and 1.95 Mo (for 3-0 
limits). 
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Fig. B.5: The obtained RVs from PARAS and TRES are plotted with re- 
spect to time. The best fit RV model with EXOFASTv2 (see Section[3.3) 
is represented by the red line, and residuals between the best fit model 
and the data are shown in the bottom panel. 
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TIC 437856897 - Sector 43 TIC 437856897 - Sector 44 TIC 437856897 - Sector 45 
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Fig. B.6: The figure represents the target pixel file for TOI-4603 in sector 43, 44 and 45 generated with tpfplotter (Aller et al.12020). The 
squared region is the aperture mask used in the photometry, whereas the size of the individual dot is the magnitude contrast (Am) from TOI-4603. 


The position of TOI-4603 is marked with ‘1’. 
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Appendix C: The corner plot showing the covariances for all the fitted parameters for the TOI-4603 
global-fit 
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Appendix D: Posterior distribution inferred for the interior modeling of TOI-4603 b 


Lumi = 2357.07*196618 
l 
Li 


Epsilon = 2.50+175 


Epsilon 


Zhomo = 0.13*002 
1 


Zhomo 


Teq = 1673.40*2273 


Teq 


Mp = 12.5892! 


Mp 


Rp = 1.057993 


Lumi Epsilon Zhomo Teq Mp Rp 


Article number, page 15 of 15 


